Two-parameter deformations of logarithm, exponential, and entropy: A consistent 
framework for generalized statistical mechanics 
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A consistent generalization of statistical mechanics is obtained by applying the maximum entropy 
principle to a trace-form entropy and by requiring that physically motivated mathematical prop- 
erties are preserved. The emerging differential-functional equation yields a two-parameter class of 
generalized logarithms, from which entropies and power-law distributions follow: these distributions 
could be relevant in many anomalous systems. Within the specified range of parameters, these en- 
tropies possess positivity, continuity, symmetry, expansibility, decisivity, maximality, concavity, and 
are Lesche stable. The Boltzmann-Shannon entropy and some one parameter generalized entropies 
already known belong to this class. These entropies and their distribution functions are compared, 
and the corresponding deformed algebras are discussed. 
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I. INTRODUCTION 

In recent years the study of an increasing number of 
natural phenomena that appear to deviate from standard 
statistical distributions has kindled interest in alternative 
formulations of statistical mechanics. 

Among the large class of phenomena which show 
asymptotic power-law behaviors, we recall anomalous dif- 
fusion [1, 2], turbulence [3, 4], transverse-momentum dis- 
tribution of hadron jets in e + e~ collisions [5], thermal- 
ization of heavy quarks in collisional process [6], astro- 
physics with long-range interaction [7] and others [8, 9]. 
Typically, anomalous systems have multifractal or hier- 
archical structure, long-time memory, and long-range in- 
teraction [10, 11]. 

The success of the Boltzmann-Gibbs (BG) theory has 
suggested that new formulations of statistical mechanics 
should preserve most of the mathematical and epistcmo- 
logical structure of the classical theory, while reproduc- 
ing the emerging phenomenology of anomalous systems. 
To this end, new entropic forms have been introduced, 
which would generalized the classical one introduced by 
Boltzmann and Gibbs and, successively, by Shannon in 
a different context (BGS entropy) 
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(1.1) 



There is no systematic way of deriving the "right" en- 
tropy for a given dynamical system. Among the many 



generalizations of the BGS entropy, one can find the en- 
tropies considered by Renyi [12], by Tsallis (q entropy) 
[13], by Abe [14], by Tsallis, Mendes and Plastino (escort 
entropy) [15], by Landsberg and Vedral [16], and recently 
by Kaniadakis (k entropy) [17, 18]. For a historical out- 
line see Ref. [19]. 

Generalizations lead to abandoning part of the original 
mathematical structure and properties. For instance, it 
is known that the BGS entropy is of trace form [20, 21] 



A' 



(1.2) 
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with A(pi) — ]n{pi); on the contrary the Renyi entropy, 
the Landsberg- Vedral entropy, and the escort entropy are 
not of trace form. The Renyi entropy and the Landsberg- 
Vedral entropy are concave only for < q < 1, while the 
escort entropy is concave only for q > 1 [22]. 

A fundamental test for a statistical functional 0(p) of 
the probability distribution to be physically meaningful 
is given by the Lesche stability condition [23]: the rel- 
ative variation [G(p) — 0(p')]/snp[0(p)] should go to 

zero in the limit that the probabilities p i — ► p[. This 
stability condition for a functional is a necessary but not 
sufficient condition for the existence of an associated ob- 
servable. Lesche [24] showed that, adopting the measure 
\\p — p'\\i = J2i \Pi ~p[ \ as estimator of the closure of the 
two distributions, the BGS entropy is stable, while the 
Renyi entropy is unstable (except for the limiting case 
q = 1 corresponding to the BGS entropy) [25-27]. Also 
the Landsberg- Vedral entropy and the escort entropy do 
not satisfy the Lesche criterion [22]. On the other hand it 
is already known that the Abe entropy [23] , the q entropy 
[28], and the k entropy [29] are stable. 

In the present paper, a natural continuation of the 
work in Ref. [20], we consider the trace- form entropy 
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given by Eq. (1.2), where A(x) is an arbitrary analytic 
function that represents a generalized version of the loga- 
rithm, while its inverse function is the corresponding gen- 
eralized exponential [30-33]. A consistent framework is 
maintained with the use of the maximum entropy (Max- 
Ent) principle. This approach yields a two-parameter 
class of nonstandard entropies introduced a quarter of a 
century ago by Mittal [34] and Sharma, Taneja [35] and 
successively studied by Borges and Roditi in Ref. [36] . It 
automatically unifies the entropies introduced by Tsallis 
[13], Abe [14], and Kaniadakis [17, 18]. 

The paper is organized as follows. In Sec. II the 
differential-functional equation for the deformed loga- 
rithm A(x) proposed in [18] is briefly reconsidered within 
the canonical ensemble formalism. In Sec. Ill we solve 
this equation obtaining the more general explicit form 
of A(x). The properties of A(x), and the consequent 
constraints on the allowed range for the deformation pa- 
rameters, are discussed in Sec. IV. The deformed algebra 
arising from the deformed logarithms and exponential is 
discussed in Sec. V. Section VI is reserved to studying 
specific members of this class: the entropies considered 
by Tsallis, by Abe, and by Kaniadakis. Other up-to- 
now overlooked cases are also discussed. The generalized 
entropies and distributions related to the deformed loga- 
rithms are studied in Sec. VII. In Sec. VIII we show that 
this two-parametric entropy is stable according to Lesche. 
We summarize the results in the final section IX. 



II. CANONICAL FORMALISM 

Guided by the form of the BGS entropy Eq. (1.1), we 
consider the following class of trace-form entropies 

N 

S(P) = -£> A (P*)' (2- 1 ) 

i=l 

with p = {pi}i=i, n a discrete probability distribution; 
one may think of A(x) as a generalization of the loga- 
rithm. 

We introduce the entropic functional 

T[ P ] = S( P )-f3' (E^- 1 ) -P (j^ EiPi ~ U ) (2 - 2) 

with 0' and Lagrange multipliers. Imposing that T\p] 
be stationary for variations of 0' and yields 

N N 
i=l i=l 

which fix the normalization and mean energy for the 
canonical ensemble. In addition, if T[p] in Eq. (2.2) 
is stationary for variations of the probabilities p j , 

/- m = 0, (2.4) 



one finds 

^ [ft A(p,)] =-/?(£,-/*), (2.5) 

where /i = —0'/ 0- 

Without loss of generality, we can express the proba- 
bility distribution p. as 

P,=«*(-f (^ -M)), (2-6) 

where a and A are two arbitrary, real and positive con- 
stants, and £(x) a still unspecified invertible function; 
we have in mind that £{x) be a generalization of, and in 
some limit reduce to, the exponential function. 

Inverting Eq. (2.6) and plugging it into Eq. (2.5), one 
finds 

_i_ MW]=A £-.(ik). (2 . 7) 

Up to this point, A(x) and £{x) are two unrelated func- 
tions and our only assumption has been that the en- 
tropy has trace form. Now if we require, by analogy 
with the relation between the exponential and the log- 
arithm functions, that £{x) be the inverse function of 
A(x), £(A(xj) = A(£(x)) = x, we obtain the following 
differential-functional equation for A(x) 

previously introduced in [18]. A simple and important 
example in this class of equations is obtained with the 
choice A = 1 and a = e _1 . In this case it is trivial to 
verify that the solution of Eq. (2.8) that satisfies the 
boundary conditions A(l) = and d A(x)/dx\ x = 1 is 
A(pj) = In pj and the entropy Eq. (2.1) reduces to the 
BGS entropy (1.1). 

In this paper we will study the deformed logarithms 
A(x) that are solutions of Eq. (2.8), the corresponding 
inverse functions (deformed exponentials), and the en- 
tropies that can be expressed using these deformed loga- 
rithms through Eq. (2.1). 

A. A counterexample 

Since Eq. (2.8) imposes a strict condition on the form 
of the function A(x), it is natural to ask what happens 
if this condition is relaxed and more general forms of 
deformed logarithms are considered. It should be clear 
from the derivation of Eq. (2.8) that, if such more gen- 
eral logarithms, which do not satisfy Eq. (2.8), are used 
to define the entropy by means of Eq. (2.1), the cor- 
responding distributions cannot be written as Eq. (2.6) 
with the "exponential" £(x) the inverse of A(x). Alter- 
natively, if one wants that the distribution be of the form 



3 



in Eq. (2.6) with the "exponential" £(x) the inverse of 
A(x), the entropy cannot be Eq. (2.1). 

For instance, let us consider the following family of 
generalized logarithms 



(2.9) 



where ln{ K }(a;) is the k logarithm [17], which we discuss 
in Sec. VI- C. This family depends on two real parameters 
k e (-1, 1) and £ > 0. We observe that Eq. (2.9) is a 
solution of Eq. (2.8), for suitable constants a and A, only 
for the case £ = 1, with hi( K ^(x) = hi{ K y(x). 
The inverse function of Eq. (2.9) is 

exp (Ki£) (a;) = exp w (sgn(a;)|2;| 1/€ ) . (2.10) 

For £ = 1 Eq. (2.10) reduces to the k exponential [17] 
exp (a;), while it reduces to the stretched exponential 
for k — 0. Moreover, this family of logarithms and ex- 
ponentials inherits from the k logarithm and the k ex- 
ponential the properties ln (jc (l/x) = — ln (K and 

Introducing the entropy 

Pi 1 

1=1 

(2.11) 

the variational principle yields 



= a exp 



(2.12) 



Equation (2.12) becomes the K-distribution [17] for £ = 
1 and the stretched exponential distribution for k = 0; 
correspondingly Eq. (2.11) reduces to the k entropy 



N 



Snip) = -J2P t m {K}O0 



(2.13) 



in the £ — > 1 limit [18], whereas in the k — > limit (a = 
A = 1) it reduces to the stretched exponential entropy 
[37, 38]: 



SdP) = E + Z> - ln ^) - + > ( 2 - 14 ) 

i=l 

where r(/i, x) is the incomplete gamma function of the 
second kind and T(fi) = r(/i, 0) is the gamma func- 
tion [39]. Equation (2.14) demonstrates that the entropy 
(2.11) has in general a form different from Eq. (2.1): 



B. Integrals of the functions A(x) and £{x) 

The fact that the generalized logarithm A(x) is a so- 
lution of the differential equation (2.8) is sufficient to 



calculate its integral 



^2 

J A{x) = x 2 A(ax 2 )-x 1 A(ax 1 ) ^ 



From the definition of the generalized exponential £(x) 
as the inverse of the generalized logarithm A(x) it is also 
simple to calculate the integral of £(x) with the change 
of variable x = £~ 1 (s) = A(s) 



£(x) dx = x 2 £{x 2 ) — X\ £(x\) 



£(x 2 )A(a£(x 2 j) -£{x 1 )A(a£(x 1 )j 



. (2.16) 



III. SOLUTIONS OF THE 
DIFFERENTIAL-FUNCTIONAL EQUATION 

In this Section we study the solutions of Eq. (2.8), 
which we rewrite in the following form: 

4-[xA(x)] -XA(-) =0 . (3.1) 
dx \aJ 

We shall select the solutions of Eq. (3.1) that satisfy ap- 
propriate boundary conditions and that keep those prop- 
erties of the standard logarithms that we judge important 
even for a generalized logarithm. 

By performing the change of variable 



x = exp 



t 

A a 



and introducing the function 



1 



A(x) = - f(Xa Inx) , 

x 



(3.2) 



(3.3) 



the homogeneous differential-functional equation of the 
first order shown in Eq. (3.1) becomes 



df(t) 
dt 



-f(t-t ) = o, 



(3.4) 



with t = Aq lna. The most general solution of 
Eq. (3.4), a differential-difference equation belonging to 
the class of delay equations [40], can be written in the 
form 

n rrii — 1 

/W = EE 0«(* 1 ,...,O* , 'e 8ii > (3-5) 
*=i j=o 

where n is the number of independent solutions s i of the 
characteristic equation 



s t - e"*° s * = , 



(3.6) 
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nii their multiplicity [s, is solution not only of Eq. (3.6), 
but also of its first m; — 1 derivatives] and multiplica- 
tive coefficients that depend on the parameters s i . In 
terms of the original function and variable the general 
solution and the characteristic equation are 



A a s. — 1 



A{x) = ^ <Hj(si, ■ ■ ■ , O Aaln(a:) 

n mi — 1 
*=i j=0 



f + K. 



A a 



(3.8) 



where m = Xasi — 1 and a'^ = (Aa) j a^ . 

In the present work, we are interested in nonoscillatory 
solutions for Eq. (3.4): this kind of solutions maintain a 
closer relation with the standard logarithm. Therefore, 
we consider only real solutions of Eq. (3.6). There exist 
four different cases depending on the value of t : (a) for 
t > we have one solution, n = 1 and m = 1; (b) for 
— 1/e < t < 0, we have two nondegenerate solutions, 
n = 2 and m i = 1; (c) for t = — 1/e we have two 
degenerate solutions, n = 1 and to = 2; (d) for t < — I/e 
there exist no solutions. 

We discuss in order the three cases (a), (b), and (c) 
that yield solutions of the delay equation (3.4) and, there- 
fore, of the corresponding Eq. (3.1). 

The case (a) is the least interesting: 



A (a;) =ai K 



(3.9) 



is just a single power and cannot change sign as one would 
require from a logarithm. 

In the case (b), we obtain a binomial solution: 

A(x) = A^k^ k 2 )x k i + A 2 (k 1 , k 2 )x k z , (3.10) 

where A 1 and A 2 are the integration constants. 

The characteristic equation (3.8) can be solved for the 
two constants a and A, 



1 



1 + S 



l/(« 1 -re 2 ) 



(3.11) 



The first condition implies A\ = —A2 = A while from 
the second one has A = 1/(k 1 — n 2 ). Equation (3.10) 
assumes the final expression 



A(x) 



X 1 — X 2 



(3.15) 



a two-parameter function which reduces to the standard 
logarithm in the [k 1 , k 2 ) — > (0, 0) limit. 

After introducing the notation A(x) = ln {jc r} (x) and 
the two auxiliary parameters k = (k 1 — k 2 )/2 and r = 
(k 1 + « 2 )/2, Eq. (3.15) becomes 



2k 

= x r \n {K} (x) . (3.16) 
The constants a and A expressed in terms of k and r 

l/(2«) 

• I I , (3.17) 



1 + r + k 



A = 



(1 + r-n) 
(1 + r + K) 



(r-K)/(2 K ) ' 



(3.18) 



are symmetric for k <-> — /c and satisfy the useful rela- 
tions (1 + r ± k) a r±K = A, which is the characteristic 
equation (3.8), and 1/A = ln {re rj (l/a), which is the dif- 
ferential equation (3.1) at x = 1. In the following, we 
call the solution (3.16) the deformed logarithm or (k, r) 
logarithm. 

Finally, we consider the case (c) . This case can be also 
obtained as the limit of case (b) when the two distinct 
solutions Sj and s 2 become degenerate. When t = —1/e, 
s = e verifies not only Eq. (3.6), but also its first deriva- 
tive 1+to exp(— tos) = 0: this solution is twice degenerate 
(si = s 2 = e); Eq. (3.7) with the boundary conditions 
(3.13) and (3.14) becomes 



A(x) = x r Inx , 



(3.19) 



where the parameter r = A a e — 1 . The standard log- 
arithm A(x) = Inx is recovered for r = 0; the same 
standard logarithm is actually recovered from the (k, r) 
logarithm in the limit (k, r) — > (0, 0) independently of 
the direction. 



(i + ^r i/( " 



-IS,) 



(3.12) 



The two arbitrary coefficients A 1 and A 2 correspond to 
the freedom of scaling x and A (a;) in Eq. (3.1). We fix 
these integration constants using the two boundary con- 
ditions 



A(1) = 
dA(x) 



dx 



= 1 



(3.13) 
(3.14) 



x=l 



IV. PROPERTIES OF DEFORMED FUNCTIONS 

The properties of the entropy (2.1), and of the cor- 
responding distribution (2.6), follow from the properties 
of the deformed logarithm A(x) = ln {K rj (x), which is 
used in its definition. Naudts [30] gives a list of general 
properties that a deformed logarithm must satisfy in or- 
der that the ensuing entropy and distribution function be 
physical. In this section we determine the region of pa- 
rameter space (k, r) where the logarithm (3.16) satisfies 
these properties and list the corresponding properties of 
its inverse, the (k, r) exponential. 
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A. (k, r) -deformed logarithm 

The following properties for the (k, r) logarithm hold 
when k and r satisfy the corresponding limitations: 



b W} (z)eC°°(iR+) ) 
d 2 , , x „ 1 



(4.1) 
(4.2) 

, (4-3) 
(4.4) 



^2 ^.rjt^ . -I«l 2 ~ 

1 

Jln {K r} (x)dx = - (i + ^ 2 _ K2 , l + r>|«| ,(4. 



dx 



(1 - r) 2 - k 2 



5) 



1 — r > |k| .(4.6) 



Equation (4.1) states that the (k, r) logarithm is an an- 
alytical function for all x > and for all k, r G K; Eq. 
(4.2) that it is a strictly increasing function for — \k\ < 
r < \k\; Eq. (4.3) that it is concave for — |k| < r < \k\, 
when | ac| < 1/2, and for —\k\ < r < 1 — when 
\k\ > 1/2; Eq. (4.4) states that the (k, r) logarithm sat- 
isfies the boundary condition (3.13); Eqs. (4.5) and (4.6) 
that it has at most intcgrable divergences for x —> + and 
x — ► +oo. These two last conditions (4.5) and (4.6) as- 
sure the normalization of the canonical ensembles distri- 
bution arising from entropy (2.1). Conditions (4.2)-(4.6) 
select the following region: 



E 2 D TZ = 



— \k\ < r < \k\ 

- 1 < r < 1 - 



if < \k\ < 



if i < k| < 1 



2 (4-7) 



which is shown in Fig. 1: every point in 1Z selects a de- 
formed logarithm that satisfies all properties (4.1)-(4.6). 
This region 1Z include values of the parameters for which 
the logarithm is finite in the limit x — > or x — > +oo. 
Note that points in region (4.7) always satisfy \n\ < 1, 
the condition obtained in Ref. [18] for the case r — 0. In 
addition, n —> implies r —> 0. Following Ref. [30] we 
introduce the dual logarithm 



In 



ln, 



(4.8) 



which is related to the original (k, r) logarithm by 



In 



*{ K ,r}( x ) =ln {«,-r}(^) = x 2r ^{ K ,r}(x)- 



(4.9) 



Note that this last result implies that Eq. (4.6) is equiv- 
alent to Eq. (4.5) on exchanging r <-> — r. Let us define 
that a logarithm is self-dual when 



In. . (x) = In* , (a;) 



(4.10) 



Then Eq. (4.9) shows that ln {re r} (x) is self-dual if and 
only if r = 0: this case coincides with the k logarithm 
[17, 18]. 



Tsallis - logarithm 
Abe - logarithm 
k- logarithm 
y- logarithm 




r 



kl 



FIG. 1: Parameter space (k, r) for the logarithm (3.16). The 
shaded region represents the constraints of Eq. (4.7) on the 
parameters. The four lines, dashed, dotted, solid, and dash- 
dotted, correspond to the Tsallis (6.1), Abe (6.5), k (6.7), and 
7 (6.12) logarithm, respectively. 



The asymptotic behavior of ln {re r} (x) for x approach- 
ing zero is 



ln {fc ,, } <» 



1 



1 



2 |k| x\ K \- r ' 



(4.11) 



in particular it results that ln {K r} (0 + ) = — oo for r < 
while, if r = \k\, ln {jc M} (0 + ) = — 1/2 is finite. Con- 
sequently, its inverse exp { ^ lKli ( x ) g° es to zero at finite 
x: the distribution function has a cut-off at finite energy, 

ex P { „,M } (-V2N) = o. 

Analogously, the behavior of ln {re r} (x) for large values of 



x is 



x \^\+r 
^+oo 2 \k\ 



(4.12) 



which implies ln {K r} (+oo) = +oo for r > — while 
again, if r = — |k|, then ln {K _| K | } (+°°) = 1/2 and 

ex P { .,-i.|}( 1 / 2 l K D = +0 °- 

The generalized logarithm verifies the following scaling 

law 



aln {K ,r}( X ) = In W a,r/„>0 B °) > 



(4.13) 



of which Eq. (4.9) is a particular case for a = — 1, and 
that becomes a\nx ~ \n(x a ) when (k, r) — > (0, 0). 

In the following we gives some useful relations. First 
the relation 

liW(zy) - \ { xr+K + xr - K ) ln {K , r} (2/) 

+ \ {y r+K + V r - K ) ^ r} (x) (4.14) 

is easily proved taking into account the definition (3.16). 
By using the identity y r ~ K = y r+K — 2k \u {k r} (y), Eq. 
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(4.14) becames 

- 2 K ln {K , r} (x) ln {K r} (y) . (4.15) 



Moreover, using 



- {x r+ «+x r -«) = -(1 + r) ln { _ } (x) + A ln {(t , p} (x/a) , 

(4.16) 

Eq. (4.14) can be rewritten as 

ln { r„ ri ( x y) = -2(l + r) ln {Kr} (x) ln {K r} (j/) 

+ ^V.Q^.d/)- (4-17) 

B. (k, r)-deformed exponential 

The deformed logarithm is a strictly increasing func- 
tion for — \k\ < r < \k\; therefore, it can be inverted for 
(k, r) £ 1Z. We call its inverse the deformed exponen- 
tial exp {K r y{ x )> whose analytical properties follow from 
those of the deformed logarithm: 



exp^^eC^J) 
d 

—— exp r , (x) > , 

A_exp {Kr} (x)>0, 
ex P { ^°) = 1 ' 



(4.18) 
(4.19) 

(4.20) 
(4.21) 



U 

/ 

— oo 


/ 



exp {K r} (x) dx = {1 + r)2 _ K2 , (4-22) 
dx 1 



exp { ^ ri (—%) (1 — r ) 2 — k2 



(4.23) 



Equation (4.18) states that the deformed exponential 
exp {K r} (x) is a continuous function for all iel, where 

I = M+, when -|k| < r < \k\, I = (-1/2 |k|, oo), 
when r = \k\, and I = (— oo, 1/2 |k|), when r = — \n\. 
Equations (4.19)-(4.23) state that exp { ^ r} (x) is a strictly 
increasing and convex function, normalized according to 
Eq. (4.21), and which goes to zero fast enough to be in- 
tegrable for x — > ±oo. 

Introducing the dual of the exponential function 



Eq. (4.9) implies 



(4.24) 



(4.25) 



which means 



(4.26) 



Only when r = does this relation reproduce that of the 
standard exponential [18]. 

The asymptotic behaviors (4.11) and (4.12) of 
In {BiP} (x) imply 



re— i-±00 



(4.27) 



for r < \k\ , 
for r > —\k\ 



ex P { „,. } ( x ) 
in particular 

eX P { K,r } (-°°) = 0+ ' 
eX P{K,r}( + 00 ) = +°° 

while 

eX P{K,r}( _1 / 2 l K D = 0+ ' When r 

exp {Kir} (+l/2|K|) = +oo , when r = -\k\ . (4.31) 
Finally, the scaling law 



(4.28) 
(4.29) 

(4.30) 



ex P {ft , r} ( x ) 



exp 



{«/a,r/, 



} {ax) 



(4.32) 



reduces to Eq. (4.26) for a = — 1 and reproduces the 
property [exp(x)] a = exp(ax) in the (k, r) — > (0, 0) 
limit. 



V. DEFORMED ALGEBRA 

Using the definition of the deformed logarithm and its 
inverse function, we can introduce two composition laws, 
the deformed sum x © y and product x ® y. 

Let us define the deformed sum: 

x@y = \n iK r} (exp {Kr} (x) exp {<i _ r} (y)) , (5.1) 
which reduces, in the (k, r) —> (0, 0) limit, to the ordi- 

0,0 

nary sum x © y = x + y. Its definition implies that the 
deformed sum satisfies the following properties: 
(a) is associative; (b) is commutative; (c) its neutral ele- 
ment is 0; (d) the opposite of x is ln {f . r) (l/ exp { ^ r} (xfj . 
If x and y are positive Eq. (5.1) yields 

h<„. ,> 0* y) = m ( », ,> (*) © ln {», r} (y) . (5-2) 

which, when (k, r) — + (0, 0), reduces to the well-known 
property log (x y) = log x + log y. 

In the same way, let us introduce the deformed product 
between positive x and y: 

x® y = exp {K _ r} (ln {Kj r} (x) + ln {re r} (y)) , (5.3) 

which reduces, for (k, r) — ► (0, 0), to the ordinary prod- 
uct x ® y = xy. This product verifies the following 
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properties: (a) is associative; (b) is commutative; (c) 
its neutral clement is 1; (d) the inverse element of x is 
exp {K „ } (-In {Bir} (x)). 

According to Eq. (5.3) we have 

CX P{„,, } ( x + y) = CX P {K ,. } ( x ) ® ex P {K ,. } (f ) ' ( 5 - 4 ) 

which reproduces in the («, r) — > (0, 0) limit the 
well-know property of the exponential exp(x) exp(y) = 
exp(x + y). 

Note that the algebraic structures A 1 = (M, ffi ) and 
A 2 = (M + , ® ) are two Abelian groups. The deformed 

sum (5.1) and product (5.3) are not distributive and the 

structure A 3 = (J? + , ffi , <S> ) is not an Abelian field. In 

any case, following Ref. [18] it is possible to define a de- 

K,r 

formed product <8> and sum © which are distributive with 

respect to © and ffi, respectively, so that the structures 

Ai = ffi, © ) and A 2 = (M + , ©, © ) are Abelian. 

Finally, from Eqs. (5.2) and (4.17), we obtain 



x ffi y = x 



ex P {K ,.}(y) 



r-\-K 



+ y 



exp r 



— 2 Kxy. 

From the practical point of view this last expression, like 
all the expressions involving the (k, r) exponentials, are 
more useful for those particular values of r and n for 
which an explicit closed form of the (k, r) exponential can 
be given. In the next section we shall see some examples. 



VI. EXAMPLES OF ONE-PARAMETER 
DEFORMED LOGARITHMS 



The relation (4.9) reads 

In (x) = — In, 

q \ / 2 — q 

while Eq. (4.15) becomes 



(6.3) 



ln(xy)=ln(x) + ln(y) + (l-q) ln» \n q {y) . (6.4) 

The g-deformed algebra already discussed in Refs. [41, 42] 
results as a particular case with r = — k=(1 — q)/2 oi 
the deformed algebra discussed in Sec. V. 



B. Abe logarithm 

As a second example we consider the constraint (r + 
l) 2 = 1 + k 2 and define q A = r + k + 1. Then the two- 
parameter logarithm in Eq. (3.16) becomes the logarithm 
associated with the entropy introduced by Abe [14] 



ln„ A (x) 



(6.5) 



which reduces to the standard logarithm for q A — ► 1. The 
invariance of Eq. (3.16) for k — > — k results in \n qA (x) 
being invariant for q A — ► l/q A . In this case the inverse 
function of the Abe logarithm (6.5), which exists because 
ln«j A (x) is monotonic for 1/2 < q A < 2, cannot be express 
in terms of elementary functions, since Eq. (6.5) is not 
invertible algebraically. We remark that Eq. (4.15) in the 
present case reads 

ln ?A (a:y) = x^- 1 \n qA (y) + yiA- 1 \n qA (x) 

+ (H 1 -^) ln 9A (z) m 9 » • (6.6) 



The two-parameter class of deformed logarithms (3.16) 
includes an infinity of one-parameter deformed loga- 
rithms that can be specified by selecting a relation be- 
tween k and r. In this section we discuss a few specific 
one-parameter logarithms that are already known in the 
literature and have been used to define entropies in the 
context of generalizations of statistical mechanics and 
thermodynamics: we show that they are in fact mem- 
bers of the same two-parameter class; we also introduce 
a few different examples of one-parameter logarithms. 



A. Tsallis logarithm 

The first example is obtained with the choice r = — k 
for —1/2 < k < 1/2. After introducing the parameter 
q = 1 + 2k [0 < q < 2] we obtain the Tsallis logarithm 
ln«j(x) = hi{( 9 _i)/2, (i- q )/2}i x ) an d the Tsallis exponen- 
tial exp q (x) as follows: 



x x ~ q - 1 
exp q (x) = [l + (l-q)x} 1 ^ 



(6.1) 
(6.2) 



C. k logarithm 

Our third example is obtained with the constraint r = 
0. Introducing the notation ln {re} (x) = \ti {k0} (x) from 
Eq. (3.16) we obtain the n logarithm and consequently 
its inverse function, namely, the n exponential introduced 
in [17, 18]: 



m M ( x ) = 



.x" — x 



exp M (x) = (^\J 1 + K 2 X 2 + KX^j 



1/re 



(6.7) 
(6.8) 



with k G (—1, 1). We remind the reader that, because 
of property (4.9), the k logarithm is the only member of 
the family that is self-dual 



(6.9) 



The function exp M (x) increases at the same rate that 
the function exp ~(—x) decreases, 



ex P {K} ( x ) ex P {fc} ( _a; ) = 1 



(6.10) 
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The K-deformed sum is obtained from the more general 
Eq. (5.1) by setting r = 0: 



x (By = x 



\/l + k 2 y 2 + y \f \ + k 2 



(6.11) 



which reduces to the ordinary sum for k —* 0. The op- 
posite approach, i.e., starting from the K-deformed sum 
(6.11) to obtain the k logarithm and k exponential has 
been taken in Refs. [17, 18], where it is shown that the 
K-deformed sum is the additivity law of the relativistic 
momenta. 



D. Other examples 

If we define the parameter w — r/\n\, we observe that 
when w = 0, ±1/3, ±1/2, ±1, ±5/3, ±2, ±3, ±5, and 
±7, the inverse function of the deformed logarithm can 
be found by solving an algebraic equation of degree not 
larger than 4; the corresponding deformed exponential 
can be written explicitly. In particular, the cases w = 
and ±1 correspond, respectively, to the k logarithm and 
to the q logarithm; the remaining cases arc different. 
Among these additional logarithms, and corresponding 
exponentials, only the cases w = ±1/3 and ±1/2 satisfy 
all the requirements discussed in Sec. IV. 

We consider explicitly the case r = ±|k|/3. Introduc- 
ing the parameter 7 = ±2 |k|/3, Eq. (3.16) defines a 
generalized logarithm 



log^ {x) = 



x 2 ~i -x- 



3 7 



(6.12) 



which reduces to the standard ones in the 7 limit. 
This logarithm is an analytical, concave, and increasing 
function for all x > 0, when —1/2 < 7 < 1/2. 
If 7 is positive the asymptotic behaviors for x 

x — > 00 are 

r 2 -y 

~ - -5— . in, 0*0 



+ 0+ 



37 



x- 

++00 37 



and 



(6.13) 



Since the logarithm (6.12) satisfies the duality relation 
In, (x) = — ln_, (1/x) , the asymptotic behaviors for 
x — > + and x — > 00 in Eq. (6.13) are exchanged when 
7 < 0. 

The corresponding 7 exponential is 

( 1 \ V3 

1 + yj\ - 47 3 a; 3 \ 



exp (x) 



+ 



1 - y/1 - 4j 3 X 3 



1/3 



1/7 



(6.14) 



which is an analytic, monotonic, and convex function for 
all x E M when —1/2 < 7 < 1/2, and reduces to the 
standard exponential in the limit 7^0. The asymptotic 
power-law behaviors for 7 positive are 

exp (x) ~ (37 1 2; I )~ 1/7 , exp (x) ~ (37a;) 1727 , 

(6.15) 



where it is clear that the asymptotic behaviors for x — > 
+00 and x — > — 00 arc exchanged when 7 changes sign, 
coherently with the property exp 7 (a;) exp_ 7 (— x) = 1. 
Finally the deformed sum given by Eq. (5.1) becomes 



x(& y = — 3jxy 
+x 



1 + y/1 - 4 7 3 ^ 



+y 



1 + y/l- 47 3 X 3 



1/3 



1/3 



1 - Jl - 4 7 3 y 3 



+ 



1 - y/i - 4 7 3 X 3 



1/3- 



1/3- 



(6.16) 



and reduces to the ordinary sum for 7^0. 



VII. ENTROPIES AND DISTRIBUTIONS 

Having obtained the deformed logarithm as a solution 
of the differential equation (3.1), the corresponding gen- 
eralized entropy follows from Eq. (2.1): 



N 



N 



»=i 



Pi _ Pi 
2k 



(7.1) 

We observe that this class of entropies coincides with the 
one introduced by Mittal [34] and Sharma and Taneja [35] 
(MST) and successively derived by Borgcs and Roditi in 
Ref. [36] by using an approach based on a biparametric 
generalization of the Jackson derivative. 

Since the entropy is defined in terms of the deformed 
logarithm (3.16), the properties (4.1)-(4.6) of ln {K r} (x) 
assure that the entropy (7.1) satisfies many of the prop- 
erties satisfied by the standard BGS entropy (1.1). In 
particular, it is (a) positive definite, S K r (p) > for p £ 
[0, 1]; (b) continuous; (c) symmetric, S K r (p) = S Kjr (q) 
with p = (pi, ...,p N ) and q = (p T (i), . . . , p T (N)) where 
r is any permutation from 1 to N; (d) expansible, which 
means that S K (p) — S K r (q) for p = {pi, . . . , p^) and 
q = (pi, . . . , p N , 0, . . . , 0); (e) decisive, in the sense 
that 5 re , r (p(°)) = where = (0, . . . , 1, . . . , 0) is 
a completely ordered state; (f) maximal, which means 
that entropy reach its maximal value when the distribu- 
tion is uniform: Max[S(p)] for p = p^ u \ with p^ = 
(1/N, 1/JV); and, finally, (g) concave. Moreover, it 
will be shown in the next section that the whole family 
of entropies (7.1) satisfies the Lesche inequality. 

We observe that for a uniform distribution p^ u \ the 
entropy (7.1) assumes the expression 



(<•"-■>) 



In 



{K,r} 



(7.2) 



and only for the case r = 0, according to Eq. (4.9), does 
it become 



S„(pW)=ln M (N), 



(7.3) 
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which is the generalization of the well-known Boltzmann 
formula and gives the entropy of a nonextensive micro- 
canonical system as the deformed logarithm of the num- 
ber of accessible states of the system. If the alternative 
form 



(7.4) 



is adopted, the entropy of a uniform distribution reduces 
to Eq. (7.3) for any values of r and k. 

From a mathematical point of view, the properties of 
the entropy (7.1) follow from those of ln {K r} (x) in the 
range x G [0, 1], while the properties of the entropy (7.4) 
follow from those of ln {re r} (x) in the range x G [1, +oo). 
This justifies our study of the properties of ln {K r} (x) in 
the whole range x G [0, +oo). 

Regarding the relationship between the entropy of a 
system and the entropies of its subsystems, additivity 
and extensivity do not hold, in general. However, it is 
possible to show that any entropy belonging to the family 
(7.1) satisfies an extended version of the additive and 
extensive property [18]. 

In fact Eq. (7.1) can be written as 



^, r (p) = -<In {( ., p} (p)> , 



(7.5) 



which expresses the entropy S K (p) as the mean value of 
ln {re r} (p). Given two systems A and B, with probability 

distributions p A and p B , we can define a joint system 



A U B with distribution p 



AuB 



pf ® pf > where the 



deformed product ® is discussed in Sec. V. From Eqs. 
(7.5) and (5.2) it follows 

S.,^UB) = S„ r (A)+S„ r (B). (7.6) 

In Fig. 2 we plot four one-parameter entropies belong- 
ing to the family of the MST entropy as functions of p 
for system with two states of probabilities p and 1 — p. 
(a) is the Tsallis entropy 



N 



p — Pi 



(7.7) 



i=i 



Notice that the entropy (7.7) is expressed in terms 
of the Tsallis logarithm as S (p) — J2iPi m g(l/Pi); 
which is different from our choice (2.1). The property 
ln g (x) = — \n 2 _ q (1/x) shows that our choice S q (p) = 
- J2iPi ln ,(^) = HiPi ^2- q {l/pi) corresponds to a dif- 
ferent labeling of the entropy q — > 2 — q. 
(b) is the Abe entropy 



(c) is the k entropy 

SM = - 



N r) 9A 

z2 p 



Pi 



1a - 9a 



(7.8) 



^ P] +K 



A.-K 



i=l 



2k 



Notice that the entropy (7.9) due to the property 
ln {Kj (l/x) = — ln {K} (x) can be written in the form 

S K (P) = T,iPi = -^P* ln { K}(K) like the 

Boltzmann-Shannon entropy. 

(d) is the 7 entropy 



JL Bi+2 7 _ D l-7 



(710) 



Entropies with the same broken-curve style yield distri- 
butions with the same-power asymptotic behavior l/x v : 
v = 1, 4/3, 2, and 4 from top to bottom; the solid curve 
shows the Shannon entropy. 

The distribution that optimizes the entropy (7.1) with 
the constraints of the canonical ensembles (2.3) is, by 
construction, 
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(7.11) 



where we recall that the deformed exponential 
exp {K r} (x) is defined as the inverse of the deformed log- 
arithm, which exists since ln {K r} (x) is a monotonic func- 
tion. The parameter fj, is determined by '^2 i p i = 1. 

In the (k, r) — > (0, 0) limit, A = 1 and a = e _1 , and 
Eq. (7.11) reduces to the well-known Gibbs distribution 



Pt =Z(p)~ 1 exp(-(3E) , 



(7.12) 



where the partition function is given by Z((3) = exp(l — 
fj,) = Y^i ex P( — The quantity exp(— (3 E) is 
named the Boltzmann factor. Analogously we can call 
cxp {K r} {—(3E/X) the generalized Boltzmann factor. 

We observe that the distribution (7.11) can not be fac- 
torized as in Eq. (7.12): the normalization constraint is 
fulfilled by fixing \i in the generalized Boltzmann factor. 

Figure 3 shows the generalized Boltzmann factors cor- 
responding to the four one-parameter entropies of Fig. 2. 
Curves with the same style have the same asymptotic 
behavior. Given this constraint and the normalization, 
the main difference between the distributions is in the 
middle region which joins the linear region and 
the Zip-Pareto region (a; 3> !)■ 



VIII. LESCHE INEQUALITY 

An important issue is whether the entropies of the fam- 
ily under consideration are stable under small changes 
of the distribution [23, 24, 29, 31]: we want to demon- 
strate that, if the two distributions arc sufficiently close, 
the corresponding relative difference of entropies can be 
made as small as one wishes. To this end, we rewrite the 
entropy (2.1): 



N 



(7.9) S( P ) = -J2p.MPi) 
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p p 




FIG. 2: Four one-parameter entropies for several values of the deformed parameter as a function of p in a two-level system: 
(a) Tsallis entropy Eq. (7.7); (6) Abe entropy Eq. (7.8); (c) k entropy Eq. (7.9); and (d) 7 entropy Eq. (7.10). Broken curves 
with the same style show entropies whose corresponding distributions have the same power-law asymptotic decay x~" , v = 1, 
4/3, 2, and 4 from top to bottom; the solid curves show the Shannon entropy. 



N 



1=1 
N >[ 



dx 



N 



i=l 



N 



Pi 



»=i 

N -\K{pJ a ) 



A 



(f) 



dec 



= -^E^ A (|)-E / oA-^-i)*, (8.1) 

1-1 1_1 -AA(0+) 

where in the second equality we have used Eq. (3.1), in 
the last equality we have made the change of variables 
s = -XA(x/a), and A(0 + ) = lim x ^ + A(x). For the 
moment we use the notation A(x) and £(x), since we do 
not need the specific form of the deformed logarithm and 
exponential. Using the fact that A _1 (x) = £{x) for the 



class of entropies under scrutiny, one finds 
S(p) 



N - AA / 0+ ) 




-£ / 


Pi - 


,J -AA(p,/«) 








/ 


Pi - 


i=l \ 




-AA(0+) 




/ [1-^(P,«)1 



(8.2) 



where in the second equality it was used X^P. = 1; that 
a£(l/A) > a£(-s/X) > p i for -1 < s < -\A(pJa), 
and the definitions [x]+ = max(x,0) and 



A(p,a)=f [ ft -af(-i) 



(8.3) 
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FIG. 3: The generalized Boltzmann factors that correspond to entropies in Fig. 2. 



From now on we revert to the notation exp {f _ r} (x) and 
ln {K r} (x). We remark that the upper limit of the integral 
s m = — AA(0 + ) in Eq. (8.2) is s m = +oo for r ^ \k\ and 
s m = l/|2/t| for r = \k\; sec Eq. (4.11). 

The definition of A(p, s), Eq. (8.3), implies that [28]: 



N 



A(p, s) - A(q, a) \< ^ \ Pi - q t \ = \\p - q\\i , (8.4) 



i=l 



and, for values of s > —A ln {K r} (1/N), 

r N 

l-iVaex P{ _ } (-£) = 2( P( -aexp {)iiP} (-^)) 



_i=l 

< E* = 1 

i=l 



J + 



(8.5) 



from which it follows that 

A{p, s) - A(q, s)\<Na ex P{ ^ } (-^) . (8.6) 

From Eq. (8.2) the absolute difference of the entropies 
of two different distributions p = { Pi } i=1 N and q = 



{ < Zi}i=i, ...,n satisfies 



S K ,M-S«,M)\ = \ J[A(p,s)-A(q,s)}ds 

A(p, s) - A(q, s) 



+ 



I 

-l 
/ 



ds 



ds 



A(p, s) - A(q, s) 
A(p, s) - A(q, s) ds . (8.7) 



Choosing -A ln {K r} (1/7V) < I < s m , by using Eq (8.4) 
in the first integral and Eq. (8.6) in the second integral 
of Eq. (8.7), we obtain 



S„, r (p) ~ S(q) < jp-qji (£+1) 
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In particular Eq. (8.8) holds for that value I that mini- 
mizes the right-hand side of Eq. (8.8), 



-A In, 



b-«l 



.9) 



{r - T \ aN 

as long as I > —A ln {re ,(1/N), which is true when 

\p-q\i<a, (8.10) 

i.e., for sufficiently close distributions, according to the 
metric ||...||i. Introducing Eqs. (8.10) and (8.9) in 
Eq. (8.8) and performing the integration using the re- 
sult (2.16), we obtain 



1 - In 



\\p - gj 

N 



and the relative difference of entropies can be written as 

s K , r (p)-s Ki M 



< F 



Ul) 

in as 

5.12) 



with 



F. 



-Q\\u N) 



IP-gh 



In 



{K,ry 



\\p - ih 



N 



because S„ 



In, 



(8.13) 



This result demonstrates that if the two distributions 
are sufficiently close the corresponding absolute differ- 
ence of entropies can be made as small as one wishes, 
since Eq. (4.5) implies that lim x _ > o+ x m {(c ri ( x ) = 0. 

In particular, the Lesche inequality for the family of 
entropies under scrutiny is valid also in the thermody- 
namic limit TV 



oo 



l T + J im f, ».r(b-?li,^) = o. 

||p-p'||->0+ JV^oo 



(8.14) 



This last result is not trivial, since the thermodynamical 
limit introduces nonanalytical behaviors that could pro- 
duce finite entropy differences between probability dis- 
tributions infinitcsimally close. We conclude this section 
by noting that Lesche stability of the (k, r) family of en- 
tropies follows also from the general proof given in [31]. 



when the ensuing distribution function is required to be 
expressed in terms of the generalized exponential through 
the natural relation (2.6). 

The solution of this equation yields the biparametric 
family of logarithms 



ln { «, r y( x ) = xT 



x n - X 



(9.1) 



the corresponding entropy [34-36] is 



JV 



l+r 



p: - p~ 

2k 



(9.2) 



This entropy is a mathematically and physically sound 
entropy when the parameters k and r belong to the re- 
gion shown in Fig. 1 and, therefore, the (k, r) logarithm 
satisfies the set of properties (4.1)-(4.6). In particular 
these entropies satisfy the Lesche stability condition. 

Distribution functions obtained by extremizing the en- 
tropy (9.2) have power-law asymptotic behaviors: such 
behaviors could be relevant for describing anomalous sys- 
tems; a comparison between several one-parameter dis- 
tribution functions is shown in Fig. 3. 

In addition, we have shown that several important one- 
parameter generalized entropies (Tsallis entropy, Abe en- 
tropy and k entropy) are specific cases of this family; 
when the deformation parameters vanish, the family col- 
lapses to the Shannon entropy. 

Our approach yielded also new one-parameter loga- 
rithms belonging to this family, whose corresponding ex- 
ponentials can be explicitly given by algebraic methods. 

There remains the question of the relevance of each 
mathematically sound entropy to specific physical situa- 
tions. In fact a wide class of deformed logarithms satisfy 
a set of reasonable mathematical properties and physical 
constraints, in particular concavity, related to thermody- 
namic stability, and the Lesche inequality, related to the 
experimental robustness. 
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